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ABSTRACT 

The dynamic response of a flexibly-mounted ring to runout of the rotating 
seat in mechanical face seal is analyzed assuming small perturbations. It is 
found that tracking ability of the stator depends only on its dynamic character- 
istics and operating conditions and is not affected by the amount of runout. 
Three different modes of dynamic response are shown and the condition for 
parallel tracking is presented. 
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NOMENCLATURE 

coefficients, see eqs. (46) and (47), i = 1,2,3, 4 
coefficients, see eqs. (32) to (34), i = 1, 2, 3 
seal clearance along center line 
equilibrium clearance 
force 


C 

C o 

F 

F 

h 

I 

I 

K 

K 

K* 

M 

M 

m 

m 

P 

R 


g 


sp 


dimensionless force, F/Sr^ 
film thickness 

stator moment of inertia about diameter 

2 3 

dimensionless moment of inertia, Ico /Sr Q 

spring constant per unit of circumferential length, K*/27rr S p 

dimensionless spring constant, K/S 

spring constant 

moment 

dimensionless moment, M/Sr Q 
stator mass 

2 2 

dimensionless mass, mu C Q /Sr o 
pressure 

dimensionless radius, r/r Q 
radial coordinate 
stator radius of gyration 
springs radial location 
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s 

t 


X,Y,Z 


x ,y ,z 

r’ J r r 


x s’ y s’ z s 


z 

Z 


o; 


y r 

y s 

*y ^ 

e 

0 

X 

V 

T 

4 > 

(p* 

$ 

CO 


2 

seal parameter, /iu>(l - R.) (r Q /C o ) 
time 

axes of inertial reference 
axes of reference fixed to rotor 
axes of reference attached to stator 
stator axial displacement 
dimensionless displacement Z/C Q 
stator angular displacement 
rotor runout 
stator nutation 

relative tilt between stator and rotor 
tilt parameter, y*r Q /C 
angular coordinate in reference 1, 2, 3 
exponent of homogeneous solution 
viscosity 

dimensionless time, cut 

angular coordinate in reference X, Y, Z 

angular location of plane of maximum film thickness 

angular location of plane z„,y o in reference X, Y, Z 

s s 

shaft angular velocity 
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Subscripts: 


1,2,3 

axes of reference 1, 2, 3 

cr 

critical 

f 

fluid 

i 

inner radius 

m 

mean radius 

o 

outer radius 

r 

rotor 

s 

stator 

sp 

support 

X,Y,Z 

inertial reference axes 

v^s 

axes of reference x ,y . 

s s 


INTRODUCTION 

Runout of the rotating seat in high speed radial face seals is almost an in- 
evitable factor which seals designers have to face. Considering the extremely 
small gap (of the order of a few micrometers) that separates the faces in a non- 
contacting face seal (fig. 1), even the slightest runout can cause clearance var- 
iation of the order of the designed separation. Moreover, since the misaligned 
seat is rotating, the effect is of a cyclic nature and can excite vibrations of the 
flexibly-mounted seal ring. 

Determination of the conditions for dynamic tracking of the misaligned 
rotating seat has been attempted by several investigators [ 1, 2, 3] . However, 
these analyses do not cover the complete range of dynamic response, and, in 
some cases, they are incomplete, treating only one degree of freedom of the 
flexibly-mounted ring. 
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The dynamic behavior of a stationary, flexibly-mounted ring facing an 
aligned rotating-seat was analyzed in a previous paper [4] . It was found that 
when the ring is slightly disturbed from its equilibrium parallel position, one 
of three modes can prevail. The initial disturbance can decay, or increase in 
time, or it can stay constant while the ring is wobbling at half the shaft fre- 
quency. The half- frequency wobble is still a matter of controversy among seal 
researchers [ 5] , and although it has been observed experimentally [6, 7] , a 
synchronous motion of the ring seems to be predominant. 

It is the objective of this paper to extend the analysis of reference [4] to 
the more general case where the rotating seat has a certain amount of axial 
runout, and to find the dynamic response of the flexibly-mounted ring under 
such conditions . 

SEAL MODEL AND COORDINATE SYSTEMS 

The seal model and various coordinate systems are shown in figure 2. The 
reference x_^, y^, z r is fixed to the rotor which is rotating at a constant angu- 
lar velocity uj about the Z axis of an inertial reference X, Y, Z. The runout of 
the rotor is represented by the constant angle y r between axes z r and Z. 
Axis y r is always pointing to the point of minimum distance between the rotor 
and plane XY. The stator has three degrees of freedom, it can move axially 
along the Z axis and rotate about the X and Y axes. The reference 
x s ,y s . z s is attached to the stator plane but is free to rotate in it, so that 
axis y s is always pointing to the point of maximum distance between the sta- 
tor and plane XY. The resulting motion of reference x ,y , z is a combina- 

s s s 

tion of nutation y g and precession ip . 

Consider a plane O g O r AB (fig. 2) at an angle <p from the plane XZ. The 
distance between points A and B, that are located radially at distance r from 
axis Z, is the local film thickness and is given by 

h = C + y g r cos(0 - ip) - y r r cos ((p - wt) (1) 
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where C is the separation O O along the seal centerline, ip is the angle be- 

S _I! ... - _ 

tween planes z g x g a °d ZX, and wt is the angle between planes z r x r and 

ZX, corresponding to the precessions of references x g ,y g ,z g and x r , y , z^, 

respectively. 

The angle (p* corresponding to the plane of maximum film thickness can 
be found from equation (1) by setting dh/d <p = 0, hence 

y g sin((p * - ip) =y r sin((p * - cot) 

or 

y sin ip - y sin cot 

tg<p* — (2) 

y cos ib - y cos cot 
's r 

The relative tilt between stator and rotor can now be defined by an angle 
y*, measured in the plane of maximum film thickness, between the axes z 

b 

and z r . Figure 3 shows the relative position of stator and rotor which can be 
described by the relative position between a new reference 1, 2, 3 and the refer- 
ence x , y , z The reference 1, 2, 3 is free to rotate about the z axis so that 

ill y 

axis 2 always points to the point of maximum film thickness and the plane 12 
is always within the plane of the stator. The film thickness h can be described 
by the relative tilt y * and the angular coordinate 6 measured from the plane 
of maximum clearance to the plane O g O r -AB. Hence, 

h = C + y*r cos 6 (3) 

where 

0 = <p - <p* (4) 

Equating (1) and (3) and using (4) we have 

y* cos (cp - <p *) =y g cos(0 - ip) -y cos (0 - wt) (5) 

Differentiating (5) with respect to <p gives 

y* sin ((p - (p*) = y sin(<£ - ip) - y_ sin (<p - wt) 
s r 


( 6 ) 
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Squaring equations (5) and (6) and adding yields 

y * 2 =y2 +y 2 - 2y g y r cos# - cot) ( 7 ) 

Equations (5) and (6) are true for any hence, 

y* cos 4>* =y g cos jp -y r cos cot (8) 

y* sin </>* =y g sin ip -y r sin cot (9) 

SYSTEMS OF FORCES AND MOMENTS 

For small angles y and y , which is the practice in seals, the equa- 
* s 

tions of motion of the stator in the inertial reference X, Y, Z are [4] . 


mZ = 2 F^ 

(10) 

Io; x = 2 M x 

(11) 

Ic^Y = 2 M y 

(12) 


where I is the stator moment of inertia about its diameter; and 2 F, 2 M 
are the proper forces and moments acting on the stator. These forces and 
moments are contributed by the fluid film between stator and rotor, and by 
the spring support on the back of the stator (fig. 1) . 

Fluid film force and moments 

The fluid film force and moments are more readily obtained in the refer- 
ence 1, 2, 3 and then transformed into the reference X, Y, Z . These force and 
moments can be found from integration of the fluid film pressure which, in 
turn, is found from a solution of the Reynolds equation 


- 8/ h 3 8E\ = 6 J w 8h +2 ah\ 
ar V 9r/ \ do at/ 


( 13 ) 


The boundary conditions of (13) are 
p = p. at r = r t 
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p = p 0 at r = r o 

The film thickness, h, is given by equation (3) from which, by using equa- 
tion (4), we find 


— = -y*r sin0 (14) 

80 m 

• • • 

— = Z + y*r cos 0 + y*0*r sin 0 (15) 

of m m 


Equations (13), (14), and (15) are based on the narrow seal approximation 
where both circumferential pressure gradient and seal curvature are ne- 
glected [8] . 

Defining a relative tilt parameter e in the form 

y *r 
- L O 


2 

and assuming small perturbations, that is, e « 1, the approach taken in [4] 
can be followed which results in the linear expressions 


F = - 2 7r/x [ - 2 - 1 r R (1 - H.) 3 a Z 
3 ^ \ q I o m ' i' o 


(17) 



where C Q is an initial separation and 
Z = Z/C Q , hence 


- - lyV (18) 

(19) 

Z is a dimensionless axial translation 
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C - + Z) (20) 

Figure 4 shows the relative position of reference 1, 2, 3 and the inertial refer- 
ence X, Y, Z. From this description and the assumption of very small tilt y 

s 

we have 



= F 


3 


( m 4 = cos <f>* - 


M„ sin 


( 21 ) 

( 22 ) 


( M f)y = sin 4>* + M 2 COS (23) 

Flexible support contribution 

In calculating the support forces only the springs are considered while 
coulomb friction of the secondary seal is neglected. Such a model properly 
describes metal bellows seals, however, its stability limit can be used as an 
upper limit for cases where friction in the secondary seal is present. This is 
because the additional friction damping increases the seals stability, provided 
the friction is not too high and causes stator hangup. 

Assuming the stator is so designed that the spring loading is balanced at 
Z = 0, y g =0. If the spring constant per unit of circumferential length is K, 
and the springs are located at r = r g p, the support reactions to a translation Z 
and tilt y s are (fig. 2) 

( F sp) z = " 2TrKr o R sp C 0 Z (24) 


and 





(25) 
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which in the interial reference is resolved into 
(“<4 ■ ( M s P ) Xs ♦ 


and 

( M sp) y - ( M sp) x sl " * 
s 


EQUATIONS OF MOTION 

Defining a seal parameter S in the form 


S = po> 


a _ R } 2 

\ C o) 


(26) 


(27) 


(28) 


where S has dimensions of force per unit area, the forces and moments can 
be normalized by 



The spring rate per unit circumferential length is normalized by 

k=£ 

s 

Summing up all the forces and moments in the inertial reference X, Y, Z by 
using equations (17) through (28), we have 
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2 = - a^y* cos (f>* - — — |y* cos cp* + - <p*J 


y* sin <p * 


- a 2 y s cos ip (30) 


_ 2a„ 

2 My = - a y* sin (p * 

cu 


y* sin (p * 


•(?•*•) 


y* cos (f> * 


- a 2 y s Sin ip (31) 


where the dimensionless coefficients a^ag.ag are 


a = — - ^ — — — ( 1 - R.) 2 
1 4 S C 1 


(32) 


a 2 = * R sp K 


(33) 




(34) 


Since the nutation angle of the stator y and its rotation angles o' 

s X 

and o>y about the X and Y axes, respectively, are very small, they can be 
treated as vectors. Hence (see fig. 2) 


a X =7 s cos * 

(35) 

Qf v = y sin rb 
-L s 

( 36 ) 

From equations (8), (9), (35), and (36) we have 

y* cos (p * = a ^ - y^ cos cut 

(37) 

y* sin <p* = o' - y sin cut 
Y r 

(38) 

and the time derivatives are 

y* cos <£>* - (p*y * sin (/>* = + coy r sin cot 

(39) 



Substituting equations (35) through (40) in equations (30) and (31); using a 

• • • 

dimensionless time t = wt for the derivatives Z, Z, a, a,y, and and 
normalizing equations (10) through (12) we finally have 


mZ + 





(41) 


Ia" x + 2a 3 o; x + ^ + a 2 )a x + a 3 a Y =Y r (a 1 cos cot - a g sin cot) 


(42) 


Ia Y + 2a 3^Y + ( a i + a 2) a Y “ = Y r ( a i sin + a 3 cos (43) 

Equations (41), (42), and (43) are the dimensionless equations of motion of 
the stator in its three degrees of freedom . When there is no runout, that is, 
y r = 0, these equations become identical to the equations of motion obtained in 
reference [ 4] as expected . Note again that the time derivatives Z, Z, ix, 
and a in equations (41) through (43) are with respect to a dimensionless time 
t where 

T = cot 

RESULTS AND DISCUSSION 

From the equations of motion (41) through (43) it is clear that for small 
perturbations the stator motion in the axial degree of freedom is uncoupled 
with the two rotations . Moreover, the solution of (41) is independent of the 
runout, y r , and is therefore identical to the solution in reference [4] for 
y r = 0. Hence, as was found in reference [4] the stator is always stable in its 
axial degree of freedom, and any small disturbance from the equilibrium sepa- 
ration C Q will decay in time. However, in order to avoid underdamping that 
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results in unnecessary oscillations of the stator, a proper combination of the 
stator mass m and spring constant K* should be selected. From equa- 
tion (41) the condition for critical damping, and hence fastest decay of any 
axial disturbance, is 



which, by equations (33), and (34), and the definition of m results in 


< mK *>cr 



(44) 


Any combination of (mK*) > (mK*) results in underdamping and should be 

L X 

avoided . Any combination (mK*) < (mK*) results in overdamping and 
may be acceptable. However, for fastest decay of axial disturbance 
mK* = (mK*) . will be the best choice. 

The homogeneous solution of equations (42) and (43) is also identical to 
that in reference [4] for the case y r = 0 . Hence, the angular stability criter- 
ion is the same as in the absence of rotor runout. The runout y r is merely a 
forcing function; thus, it affects only the particular solution of the angular equa- 
tions of motion. The particular solution for the angles and Oy can be 
sought in the form 


Q! x = A^ cos cot + A 2 
«¥ = -Ag cos cot + A 4 


sin cot 
sin cot 


(45) 

(46) 


Substituting equations (45) and (46) into equations (42 and (43) and solving sepa- 
rately for the coefficients of sin cot and cos cot gives 
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A. = A. =y 
1 4 ' r 


a 3 + < a l + a 2 ~ 1 > a l 

4 + < a l + a 2 - T » 2 


A n = -A„ =' 


a i a 3 ~ ( a l + a 2 ~ r ^ a 3 
a 3 + ( a i + a 2 -I) 2 


Thus, the particular solution for a x and a Y is 


a x = cos o ) t + A 2 sin cot 


oiy = A^ sin cot - A 2 cos cot 


From equations (49) and (50) it is clear that this solution yields 


4 + q,2 = a 2 + a 2 


(47) 


(48) 


(49) 

(50) 


and, hence, by equations (35) and (36) the contribution of the particular solu- 
tion to the stator tilt, y , is 



(51) 


Substituting equations (47) and (48) in (51) we find 


-. 1/2 




a 2 t a 2 

a l a 3 


a 3 + (a i + a 2 ' r > 2 


(52) 


Figure 5 is a graphical representation of the stator angular motion contributed 
by the particular solution of the angular equations of motion. The complete 
motion is a superposition of the particular and homogeneous solutions . The 
homogeneous solution has the form 
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a X ~ a Xo e 
a Y = or Yo e 


Xt 

Xt 


where a XQ , a Yo is an initial disturbance. In reference [4] it was found that 
three modes of response to the disturbance ojyQ are possible depending 

on the shaft speed go. If w is less than a critical speed u> , then the initial 
disturbance a xo’ a Yo decays i° time. Hence, for co < w cr the homogeneous 
contribution disappears after a while and the stator motion is described by the 
particular solution (eqs. (46) through (49) and fig. 5). If o> > co the initial 
disturbance increases in time and the seal may fail due to stator rotor touch- 
down. At the critical speed, cj = u> cr , the homogeneous solution becomes an 
oscillatory motion with a constant amplitude and at a frequency which is half 
the shaft speed. The critical speed, co , at which this half frequency wobble 
occurs was found in [4] from the condition 


I = 4(a 2 + a x ) 

which by equations (32), (33), and the definition of I results in 


3 R 2 

9 TIT A 

2 o_m 

cr T 



r v - R Ai - p o»| 


( 53 ) 


For a circular- ring stator the moment of inertia about a diameter is 



2 

g 


(54) 


where r is the stator radius of gyration. Using equation (54), equation (53) 

b 

can be rearranged in the form 



16 



= 4 + 2 7r(l - R.) 


2 <Pj - P Q ) r o Am 


fc 

r„ 


(55) 


K* — 


spy 


Figure 6 presents stability limits based on equation (55) for seals of var- 
ious ratios R i = r i /r Q . If the seal operates in the stable regime the stator 
will track the rotor in a synchronous mode (see eqs. (49 ( and (50)) and with a 

constant amplitude y „ related to the runout y by equation (52). At condi- 

b r 

tions of critical stability, half-frequency wobble is superimposed on the basic 
synchronous tracking. When the seal operates in the unstable regime tracking 
becomes impossible and failure occurs. 

The synchronous tracking in the stable operation regime is not necessarily 
in phase with the runout. From figure 5 we see that the phase angle cot - ip 
is given by 

A 1 

COS (cot ~ Ip) - — 

Yc 


Hence, by equations (47) and (52) 


cos (cot - 4>) =' 


ig + (a 1 + a 2 - I)a x 


(' 


2 2 
a i + a 3 


,\l/2 f 

7 a 3 


a 2+ (a 1 + a 2 -I) 2 


1/2 


(56) 


Equations (52) and (56) give the amplitude y of the tracking stator and its 

phase angle (cot - ip), respectively, for a given seal geometry, spring constant, 
and operation conditions. 

For minimum leakage it is important that the stator remains parallel to 

the rotor all the time. This requires y =y , ip = cot (see fig. 2). From 

s r 

equation (56) the condition for ip - cot is 
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a 2 = I (57) 

From equation (52) we see that condition (57) also fulfills the requirement 

y =y . Hence, the condition for parallel tracking is by (57) and (33), and 
s r 

the definition of I 


Io) 2 = TrrJpH 


Substituting equations (54) in 


lr sp/ K * 


(58) 

(58), the condition for parallel tracking becomes 

(59) 


The preceding analysis is based on the assumption of small perturbations, 
e 2 « 1, which enables the linearization of the equations of motion. It is also 
assumed that the fluid film in the sealing gap is complete and does not cavitate. 
These assumptions permit analytical treatment of the complex problem of seals 
dynamics. In many applications where the runout is very small and the pres- 
sure differential is sufficient to suppress cavitation the above mentioned assump- 
tions are very close to reality. Hence, the present results may be used as a 
tool for better understanding and design of non-contacting mechanical face seals. 
In particular, stability maps similar to these in figure 6 are useful with the addi- 
tional conditions for parallel tracking represented by equation (59), and the cri- 
tical damping given by equation (44) . 

It should also be noted that the present analysis deals with a flat stator that 
has very little static stability when the high pressure is on the outer diameter 
[9] . In these cases it would be benefitial to machine the stator with a certain 
amount of coning [10] to increase its static stability. The coning also helps in 
obtaining the desired equilibrium clearance C . Except for altering the sta- 
bility maps, coning will not change the basic nature of the dynamic response. 



18 


This will still depend on the regime of the seal's operation. That is, in the sta- 
ble regime synchronous tracking at some constant amplitude; in the unstable re- 
gime failure due to inability to track; and at the transition between these two 
regimes, half frequency wobble superimposed on the basic synchronous track- 
ing. Indeed all these three modes of dynamic response were observed experi- 
mentally [7] , verifying qualitatively the present theory. 

A question often brought up among seals researchers and designers is how 
much runout can be tolerated before the stator loses its capability to track? 

The results of this work show that tracking ability of the stator is not affected 
by the amount of rotor runout but depends on the dynamic characteristics of the 
stator, the pressure differential, and the shaft speed. The runout itself is 
merely the amplitude of a forcing function and hence, does not affect the sta- 
bility of the system. A word of caution should, however, be added at this 
point. Although the amount of runout does not affect the tracking ability of the 

stator it does affect its tilt y and phase angle cot - ip. Hence, the amount of 

S' 

runout affects the relative tilt y* and, therefore, the minimum film thick- 
ness between rotor and stator (see fig. 3). Thus, a situation may exist where 
the stator steadily tracks the rotor runout but the minimum film thickness be- 
comes dangerously small, resulting eventually in seal failure. From equa- 
tion (3) we see that the minimum film thickness, corresponding to r = r Q 
and 0 = 7r, is 

h . = C - y*r (60) 

Hence, when y*r Q /C approaches unity, h min approaches zero. A situation 
like this cannot be handled by the present analysis which is based on the as- 
sumption of small perturbation e « 1 and, therefore, will require a treat- 
ment of the more general nonlinear equations of motion . Some idea on the 
effect of the runout y r on the minimum film thickness can, however, be ob- 
tained from the present results. Consider equations (7), (52), and (56) we find 
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fry 1 * 

+ + a 2 - I) 

From equation (61) we see again that the condition for parallel tracking, 

7 * = 0, is met when a 2 = I. In any other case a 2 ^ I, after ensuring steady 
tracking (operation in the stable regime, see fig. 6 ) and critical damping, one 
can decide on a maximum permissible relative tilt y* and use equation (61) 
to calculate the amount of runout y r that will be damaging. 

It is also worthwhile noting that the presence of friction in stators having 
secondary seal O ring and antirotation locks may prevent parallel tracking at 
all. In such seals dynamic instability can become a lesser problem (because 
of the additional damping) compared with the possibility of excessive relative 
tilt y* caused by rotor runout. 

CONCLUDING REMARKS 

The dynamic response of a flexibly-mounted stator to runout of the rotating 
seat in radial face seals is analyzed. It is found that successful tracking is not 
affected by the amount of runout but rather by the nature of stator dynamic 
characteristics. Three different modes of dynamic response are possible. If 
the seal operates in a stable regime, synchronous tracking at constant ampl i- 
tude prevails . In the unstable operating regime tracking is impossible, re- 
sulting in face contact due to angular instability. At the transition between 
these two regimes half frequency wobble is superimposed on the basic syn- 
chronous tracking. A condition for parallel tracking was found which is 
important for minimum leakage. It is also shown that the amount of runout, 
although not affecting the dynamic response of the system, may be responsible 
under non-parallel tracking conditions to seal failure due to an excessive rela- 
tive tilt between stator and rotor. 

The analysis is based on certain simplifying assumptions but nevertheless 
provides a good insight into the dynamic behavior of non-contacting face seals. 
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Figure 1. - Schematic of a radial face seal. 
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Figure 2. - Seal model and coordinate systems. 




Figure 3. - Relative position between stator and rotor. 



Figure 4. - Orientation of inertial system 
X.YjZ in rotating system 1, 2, 3. 
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Figure 5. - Graphical representation of the 
particular solution for q x and cy 



Figure 6. - Stability maps for the seal stator at various radius ratios Rj. 
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